A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 
06(02.01.2, 02.02.1, 08.14.1, 02.08.1) 



ASTRONOMY 
AND 
ASTROPHYSICS 
2.2.2008 



O 
O 

(N 

> 
O 

o 



Pseudo-Schwarzschild description of transonic spherical 
accretion onto compact objects 

Tapas K. Das 1 and Aveek Sarkar 2 

1 Inter University Centre For Astronomy And Astrophysics, Post Bag 4 Ganeshkhind, Pune 411 007, India 
email: tapas@iucaa.ernet.in 

2 National Centre For Radio Astrophysics, TIFR, Post Bag 3 Ganeshkhind, Pune 411 007, India 
email: sarkar@ncra.tifr. res. in 

Received 6 March 2001/ Accepted 9 May 2001 



Abstract. A number of 'modified' Newtonian potentials of various forms are available in the literature which ac- 
curately approximate some general relativistic effects important for studying accretion discs around a Schwarzschild 
black hole. Such potentials may be called 'pseudo-Schwarzschild' potentials because they nicely mimic the space-time 
y—i • around a non-rotating/slowly rotating compact object. In this paper, we examine the validity of the application of 
some of these potentials to study the spherically symmetric, transonic, hydrodynamic accretion onto a Schwarzschild 
■ black hole. By comparing the values of various dynamical and thermodynamic accretion parameters obtained for flows 
\ using these potentials with full general relativistic calculations, we have shown that though the potentials discussed in 
£SJ ■ this paper were originally proposed to mimic the relativistic effects manifested in disc accretion, it is quite reasonable 
t-H ' to use most of the potentials in studying various dynamical as well as thermodynamic quantities for spherical accre- 
tion to compromise between the ease of handling of a Newtonian description of gravity and the realistic situations 
| described by complicated general relativistic calculations. Also we have shown that depending on the chosen regions 
of parameter space spanned by specific energy £ and adiabatic index 7 of the flow, one potential may have more 
Or importance than another and we could identify which potential is the best approximation for full general relativistic 
Q flow in Scwarzschild space-time for particular values of £ and 7. 

a: 

c/2 ■ Key words: accretion, accretion discs — black hole physics — hydrodynamics 
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1. Introduction 

Stationary, spherically symmetric and transonic hydrodynamic accretion of adiabatic fluid on to a gravitating astro- 
physical object at rest was studied in a seminal paper by Bondi (1952) using purely Newtonian Potential and by 
including the pressure effect of the accreting material. Later, Michel (1972) discussed fully general relativistic poly- 
tropic accretion on to a Schwarzschild black hole by formulating the governing equations for steady spherical flow of a 
perfect fluid in the Schwarzschild metric. Following Michel's relativistic generalization of Bondi's treatment, Begelman 
(1978) discussed some aspects of the critical (sonic) points of the flow for such an accretion. Using an unrelaxed mono- 
energetic particle distribution and assuming the fact that the relaxation time of such a particle distribution is very 
long compared to the typical flow time scale or dynamical time scale of steady accretion on to black holes, Blumcnthal 
and Mathews (1976) developed a model where the connection between the nonrelativistic to the relativistic regime of 
the spherically accreting material could be established. Taking the fully ionized one-temperature (T e ieutron — T pro ton) 
hydrogen gas (governed by an exact relativistic equation of state) to be the fundamental constituent of the accret- 
ing material, Brinkmann (1980) treated spherically symmetric stationary accretion in Schwarzschild space time and 
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showed that the temperature of accreting material at the Schwarzschild radius is one order of magnitude smaller than 
the flow temperature obtained by using a simple polytropic equation of state. Recently, Malec (1999) provided the 
solution for general relativistic spherical accretion with and without back reaction and showed that relativistic effects 
enhance mass accretion when back reaction is neglected. 

Meanwhile, the theory of the accretion disc found prior importance because of the fact that disc accretion describes 
more realistic astrophysical situations found in nature. The beginning of modern accretion disc physics is traditionally 
attributed to the two classical articles by Shakura & Sunyaev (1973) and Novikov & Thorne(1973) . While Shakura 
& Sunyaev (1973) calculated the disk structure and related phenomena using purely Newtonian potential, Novikov & 
Thorne provided a fully general relativistic description of accretion discs around black holes; later on, some aspects of 
which were modified by Riffert and Herold(1995). However, rigorous investigation of transonic disc structure was found 
to be extremely complicated in full general relativistic space time (Chakrabarti (1996) and references therein). At the 
same time it was understood that as relativistic effects play an important role in the regions close to the accreting 
black hole (where most of the gravitational potential energy is released), purely Newtonian gravitational potential 
(in the form $ Newton = —^-) cannot be a realistic choice to describe transonic black hole accretion in general. To 
compromise between the ease of handling of a Newtonian description of gravity and the realistic situations described 
by complicated general relativistic calculations, a series of 'modified' Newtonian potentials have been introduced to 
describe the general relativistic effects that are most important for accretion disk structure around Schwarzschild and 
Kerr black holes ( see Artemova et. al.(1996) for further discussion). Introduction of such potentials allows one to in- 
vestigate the complicated physical processes taking place in disc accretion in a semi-Newtonian framework by avoiding 
pure general relativistic calculations that most of the features of spacetime around a compact object are retained and 
some crucial properties of the analogous relativistic solutions of disc structure could be reproduced with high accu- 
racy. Hence, those potentials might be designated as 'pseudo-Kerr' or 'pscudo- Schwarzschild' potentials, depending 
on whether they are used to mimic the space time around a rapidly rotating or non rotating/ slowly rotating (Kerr 
parameter a ~ 0) black hole respectively. 

It is important to note that although a number of such 'pseudo' potentials are available in the literature to study 
various aspects of disc accretion, no such potentials are available which had been solely derived to describe spher- 
ically symmetric accretion on to a Schwarzschild (or Kerr) black hole. In this paper, we will concentrate on some 
of the 'pseudo-Schwarzschild' disc potentials (potentials introduced to study accretion discs around a Schwarzschild 
black hole) to investigate whether those potentials could be used to study Spherical accretion, and if so, how 'good' 
the choice would be for various such potentials. Also, we would like to check which potential among those would 
be the 'best-fit' to approximate the full general relativistic description of transonic, spherically symmetric accretion 
on to a Schwarzschild black hole. In doing so, we solve the equations of motion of spherically accreting fluid in full 
Schwarzschild space-time as well as for motion under various 'pseudo'-potentials, to study the variation of different 
dynamical and thermodynamic quantities (like Mach number of the flow, flow temperature etc.) with radial distance 
measured from the accreting black hole for the full general relativistic spherical flow (hereafter FGRSF) as well as for 
accretion using various 'pseudo- Schwarzschild' potentials. We then compare the results obtained using such poten- 
tials with the solutions of exact relativistic problems in a Schwarzschild metric. The plan of the paper is as follows: 
In next section, we will describe four 'pseudo-Schwarzschild' disc potentials available in the literature and some of 
their basic features. In §3, we will provide the basic equations governing spherically symmetric accretion in full rel- 
ativistic as well as in various 'pseudo'-relativistic spacetimes. In §4 we will discuss how to solve those equations to 
find various dynamical quantities which are to be mutually compared and we present our results. Finally in §5 we 
conclude by discussing the suitability of various 'pseudo' potentials in approximating the results obtained from exact 
relativistic calculations. For the rest of this paper, we will use the terms 'modified- Newtonian potential' and 'pseudo 
(Schwarzschild) potentials' synonymously. 

2. Some basic features of various 'pseudo-Schwarzschild' potentials 

From now, we will define the Schwarzschild radius r g as 

_ 2GM BH 

(where Mbh is the mass of the black hole, G is universal gravitational constant and c is velocity of light in vacuum) so 
that the marginally bound circular orbit r\> and the last stable circular orbit r s take the values 2r g and 3r g respectively 
for a typical Schwarzschild black hole. Also, total mechanical energy per unit mass on r s (sometimes called 'efficiency' 
e) may be computed as —0.057 for this case. Also, we will use a simplified geometric unit throughout this paper where 
radial distance r is be scaled in units of r g , radial dynamical velocity u and polytropic sound speed a of the flow is 
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scaled in units of c (the velocity of light in vacuum) , mass m is scaled in units of Mbh and all other derived quantities 
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Fig. 1: Newtonian potential &Newton(r) and other pseudo-potentials <&»(r)'s (i = 1,2,3,4) are plotted as a function of the 
logarithmic radial distance from the accreting black hole. Except $2, all other pseudo-potentials have a singularity at one 



Schwarzschild radius r 5 
for details. 



For r > 2r 9 , while the relative stiffness factor S is maximum for $4, it is minimum for $2, see text 



would be scaled accordingly. Also, for simplicity, we will use G = c = M = 1. Below we would like to briefly describe 
four different 'pseudo- Schwarzschild' potentials (expressed in the system of units discussed above) and to provide the 
'free fall' acceleration obtained using such potentials in compact form. 
Paczyhski and Wiita (1980) introduced a 'pseudo-schwarzschild' potential of the form 



*i = 



1 



2(r-l) 



(la) 



which accurately reproduces the positions of r s and r/, and gives the value of efficiency to be —0.0625. Also the 
Keplarian distribution of angular momentum obtained using this potential is exactly same as that obtained in pure 
Schwarzschild geometry. Although this potential (as well as the other 'pseudo' potentials available in the literature) 
does not satisfy the boundary condition exactly on the horizon, however, close to the horizon, accretion flows are 
supposed to be highly supersonic and the dynamical infall time scale becomes too small to allow significant radiation 
from infalling fluid. Even if there are significant amounts of radiation, it does not contribute too much to the overall 
radiation coming out from the disc, due to the fact that radiation emitted from the very close viscinity of the black hole 
is highly rcdshiftcd. Thus, errors in calculation close to the horizon may have very little impact on the emitted disc 
spectrum and for all practical purposes, the Paczyhski and Wiita (1980) potential is considered to be the best approx- 
imation of pure Schwarzschild spacetime while treating the disc accretion (Artemova 1996) . It is worth mentioning 
here that this potential was first introduced to study a thick accretion disc with super Eddington Luminosity. Also, it 
is interesting to note that although it had been thought of in terms of disc accretion, it is spherically symmetric (with 
a scale shift of r g ). So in principle, it can definitely be used to study spherical accretion also on to a Schwarzschild 
black hole. 

To analyse the normal modes of accoustic oscillations within a thin accretion disc around a compact object (slowly 
rotating black hole or weakly magnetized neutron star), Nowak and Wagoner (1991) approximated some of the domi- 
nant relativistic effects of the accreting black hole (slowly rotating or nonrotating) via a modified Newtonian potential 
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of the form 



1 3 f I 



) 



2 



(16) 



$2 has correct form of r s as in the Schwarzschild case but is unable to reproduce the value of r&. This potential has 
the correct general relativistic value of the angular velocity (as measured at infinity) at r s . Also it reproduces the 
radial epicyclic frequency k (for r > r s ) close to its value obtained from general relativistic calculations. However, this 
potential gives the value of efficiency as —0.064 which is larger than that produced by <j>i, hence the disc spectrum 
computed using $2 would be more luminous compared to a disc structure studied using $1. 

Remembering that the free-fall acceleration plays a very crucial role in Newtonian gravity, Artemova et. al. (1996) 
proposed two different 'pseudo' potentials to study disc accretion around a non-rotating black hole. The first potential 
proposed by them produces exactly the same value of the free-fall acceleration of a test particle at a given value of r 
as is obtained for a test particle at rest with respect to the Schwarzschild reference frame, and is given by 



The second one gives the value of the free fall acceleration that is equal to the value of the covariant component of 
the three dimensional free-fall acceleration vector of a test particle that is at rest in the Schwarzschild reference frame 
and is given by 



Efficiencies produced by $3 and $4 are —0.081 and —0.078 respectively.The magnitude of efficiency produced by $3 
being maximum, calculation of disc structure using $3 will give the maximum amount of energy dissipation and the 
corresponding spectrum would be the most luminous one. However, as both $3 and $4 stems from the consideration 
of free fall acceleration and calculates the dependence of free fall accleration on radial distance in the Schwarzschild 
metric (which describes a spherically symmetric gravitational field in vacuum), it appears to be quite justified to use 
those potentials to study spherically symmetric accretion. 

From now we will refer to all these four potentials as $i in general where {i — 1, 2, 3, 4} would correspond to $1 (eqn. 
1(a)), $2 (eqn. 1(b)), $3 (eqn. 1(c)) and $4 (eqn. 1(d)) respectively. In Fig. 1, we plot various <!>i's as a function of 
the radial distance measured from the accreting black hole in units of r g . Also in the same plot, purely Newtonian 
potential $ Newton is plotted. If we now define a quantity Si to be the 'relative stiffness' of a potential <j>j as: 



(that is, Si is a measure of the numerical value of any zth potential at a radial distance r), we find that for r > 2r g : 



which indicates that while $2 is a 'flatter' potential compared to the pure Newtonian potential <& Newton, all other 
'pseudo' potentials are 'steeper' to $ Newton for r > 2r g . 

One can write the modulus of free fall acceleration obtained from all 'pseudo' potentials except for $2 in a compact 
form as 




(lc) 




(Id) 




S2 < SNewton < Si < S3 < S4 



2 r 2-*i(r- l) 5 ' 

I and 64 = 1. <& i denotes the absolute value of the space derivative of i.e. 



(2a) 



where 6± — 



2, £ 3 - 




d^ 



dr 



whereas acceleration produced by <j> 2 can be computed as, 




(26) 



In the next section, we would like to describe how one can investigate transonic spherical accretion using these po- 
tentials. Also, we will discuss how to calculate various dynamical quantities for full general relativistic bondi flow in 
a Schwarzschild metric. One standard method to investigate classical transonic bondi flow is to formulate the basic 
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conservation equations, i.e., conservation of baryon number (obtained by integrating continuity equation) and con- 
servation of specific energy (obtained by integrating Euler's equation using a specific equation of state for accreting 
matter) and then to simultaneously solve these conservation equations to get critical (sonic) quantities as functions 
of various accretion parameters (like specific energy £, adiabatic index 7 or accretion rate M in in the flow) and also 
to calculate the values of various dynamical and thermodynamic quantities (like Mach number M, of the flow, flow 
temperature, T, etc.) as functions of various accretion parameters (like £,7, Mi„ etc.) or radial distance (measured 
from the central accretor in the unit of Schwarzschild radius r g ). Also a common practice is to study the variation of 
Mach number M (ratio of the local dynamical velocity u to the local sound velocity a; M = -) with radial distance 
r (measured in units of r g ) to investigate the 'transonicity' of the flow for a fixed set of input parameters. We will 
perform the above mentioned calculations for accretion in pure Schwarzschild spacetime as well as for accretion in 
'pseudo-schwarzschild' space time using <E>j's and compare the results obtained for various <!Vs with exact relativistic 
calculations. As accretion onto a black hole is necessarily transonic to satisfy the inner boundary condition at the 
event horizon, unlike Bondi's (1952) original work, here we will concentrate only on transonic solutions, i.e., we will 
deal with accretion (wind) which is subsonic (supersonic) far away from the black hole and approaches (moves away 
from) the hole supersonically (subsonically) after crossing a sonic point (the location of which can be determined as 
a function of £ and 7, see §3) on its way towards (away from) the accretor. 

3. Governing equations 

3.1. Full General Relativistic Spherical Flow (FGRSF) in Schwarzschild space-time 
For a schwarzschild metric of the form 

ds 2 = dt 2 (l - ^ - dr 2 (l - J) 1 - r 2 (d9 2 + sin 2 9d^ 2 ) 

the energy momentum tensor T a/3 for a perfect fluid can be written as (Shapiro & Teukolsky 1983) 

T a(3 = ea Q / + p (u a u fi - g a ?) 

where e and p are proper energy density and pressure of the fluid (evaluated in the local inertial rest frame of the 
fluid) respectively and u a is the four velocity commonly known as 

dr a 
ds 

Equations of motion which are to be solved for our purpose are, 

1) Conservation of mass flux or baryon number conservation: 

(pu a ). a =0 (3a) 

and 

2) Conservation of momentum or energy flux (the general relativistic Euler equation obtained by taking the four 
divergence of T af3 ): 

(e + p) u a . p u fJ = -p ia - u a p^u p (36) 
where the semicolons denote the covariant derivatives. 

Following Michel (1972), one can rewrite eqn. 3(a) and eqn. 3(b) for spherical accretion as 

Airpur 2 = M in (4a) 

and 



P + e \ {, 1 , . 2 



p j{ 1 - 1 r +u2 )= C (45) 

as two fundamental conservation equations for time independant hydrodynamical flow of matter on to a Schwarzschild 
black hole without back-reaction of the flow on to the metric itself. Mi n is the mass accretion rate and C is a constant 
(related to the total enthalpy influx) to be evaluated for a specific equation of state. 
For a polytropic equation of state i.e, 

p = Kp 1 



(i 
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and defining M. = Mi„'y n K n , where A4 is a measure of the entropy of the flow (n is a polytropic constant of the flow 
defined as n — (7 — 1) ) and is another conserved quantity of the flow, one can rewrite conservation equations 4(a) 
and 4(b) as (Chakrabarti (1996) and references therein): 

5 = l + ^ = 1+ (^)(^i) 2 (5a) 

and 

2 



M = 4tt uu t r z (56) 

\ 1 — na z J 

Where £, h, Ut and a are the conserved specific energy of the flow excluding its rest mass, specific enthalpy, specific 
binding energy and local adiabatic sound speed respectively. Eq. 5(b) may be considered as the outcome of the 
conservation of mass and entropy along the flow line. The expression for a can be written as (Weinberg (1972), Frank 
et. al. (1992)): 



Constant Specific Entropy 




(6) 



Where T is the flow temperature, \i is the mean molecular weight and m#^m p is the mass of the hydrogen atom. 
One can now easily derive the expression for the velocity gradient (4j:) (by differentiating eqn. 5(a) and 5(b)) as 

du u(l-u 2 ){a 2 (4r-3)-l} 

dr 2r (r - 1) (u 2 - a 2 ) 1 ° j 

Since the flow is assumed to be smooth everywhere, if the denominator of eqn. 7(a) vanishes at any radial distance 
r, the numerator must also vanish there to maintain the continuity of the flow. One therefore arrives at the so called 
'sonic point (alternately, the 'critical point' El) conditions' by simultaneously making the numerator and denominator 
of eqn. 7(a) equal zero. The sonic point conditions then can be expressed as follows 

(76) 



4r c - 3 

For a specific value of £ and 7, location of the sonic point r a J can be obtained by solving the following equation 
algebraically 



(7 - 1) (4r c - 3) / 3 
^4, c (7-l)-(37 + 4)V 1 ^ + 1 = ° (7C) 

The spherical surface of radius r = r 9 J can be defined as the 'accoustic horizon' because for r < r 9 c r , u > a and 
any accoustic disturbances created in this region are advected towards the black hole. Thus no accoustic disturbances 
created within this radius can cross the accoustic horizon and escape to the region r > r c . 

To determine the behaviour of the solution near the sonic point, one needs to evaluate the value of (%) at that point 
(we denote it by (4 s ) ) by applying L 'Hospitals' rule to eqn. 7(a). It is easy to show that (4jr) can be obtained by 
solving the following quadratic equations algebraically: 

du\ 2 (7 - 1) (I6rg - 16r e - 8 7 r e + 67 + 3) f du 



= (7d) 



dr Jc 3r c (4r c -3)^ \ dr 

(7 - 1) (2r c - 1) (24r 2 - 28r c - 8^7 + 4r c7 + 37 + 6 
2r c (4r c -3) f (4r c -3) (r c - 1) 



It is now quite straightforward to simultaneously solve eqn. 5(a) and eqn. 5(b) to get the integral curves of the flow 
for a fixed value of £ and 7. Detailed methodology for this purpose will be discussed in §4. 



Hereafter, we will use 'critical ponts' and 'sonic points' synonymously. 
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3.2. Spherical flow in various Pseudo-Schwarzschild space time 

For any 'pseudo-Schwarzschild' potential <&i, the equation of motion for spherically accreting material onto the accretor 
is given by 

du du 1 dP ^ i ,„ , 

at or p or 

where symbols have their usual meaning. The first term of eqn. 8(a) is the Eulerian time derivative of the dynamical 
velocity at a given r, the second term is the 'advective' term, the third term is the momentum deposition due to 
pressure gradient and the final term is due to the gravitational acceleration for a particular ith potential The 
continuity equation can be written as 



dp 1 d , 
dt r 2 dr 







(86) 



For a polytropic equation of state, the steady state solution (apart from a geometric factor of Air) of eqn. 8(a) and 
eqn. 8(b) is 

1) Conservation of specific energy £ of the flow: 



7-1 



(9a) 



and 



2) Conservation of Baryon number (or accretion rate M in ): 

M in = At: pur 2 

Using M as defined earlier (M = M in -f n K n ), eqn. 9(b) can be rewritten as 

M = Aita 2n ur 2 



(96) 



(9c) 



It is now quite straightforward to derive the space gradient of dynamical velocity for flow in any particular ith 

potential <J>i as 

'rlii\ cb- 

(10a) 



d:r 



u — — 



Using either eqn. 2(a) or eqn. 2(b), one can substitute the value of any specific ith in eqn. 10(a) to get the value of 
( eF ) • f° r m °tion under particular ith pseudo-Schwarzschild potential. One can also calculate the sonic point quantities 
(as described in §3.1) as 



— *i' 
2 



where superscript i denotes the specific value of sonic quantities for a particular ith potential <&j, and $ i 



(106) 



is the value 



of (^ff) evaluated at the corresponding sonic point r l c . The value of sonic point r\ for any ith potential <&i can be 
obtained by algebraically solving the following equation 



2 Vl-l 







(10c) 



where $, 



is the value of ith potential at the corresponding sonic point r*. Similarly, the value of for any $j at 



its corresponding sonic point r* can be obtained by solving the following quadratic equation: 



+ 2.829 (7-I) 



\ 



du\ 

drj ci r\ 



= 



(lOd) 
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where $; 



is the value of 4^ at the corresponding critical point r l c . 



One can simultaneously solve eqn. 9(a) and eqn. 9(b) (alternatively, eqn. 9(a) and eqn. 9(c)) for any specific $^ for 
a fixed value of £ and 7 to obtain various dynamical and thermodynamic flow quantities. We discuss the detailed 
methodology in §4. 
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Fig. 2 (a,b,c,d): The deviation (Arl) of the value of sonic points r*'s obtained for flows using various $i's from rf as a 
function of specific energy £ and adiabatic index 7 of the flow is plotted for flows in $1 (Fig. 1(a)), $2 (Fig. 1(b)), $3 (Fig. 1(c)) 
and $4 (Fig. 1(d)). It is observed that while Arl is a minimum for ultra-relativistic (7 = |) flow, Ar* becomes a minimum as 
the flow approaches a pure-non-relativistic (7 = |) limit, see text for detail. 



4. Solution procedure and results 
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4-1- Sonic points as a function of £ and 7 

Solving eqn. 7(c) and eqn. 10(c) for a fixed value of £ and 7, the sonic point could be obtained for FGRSF and for 
flows in various pseudo-potentials respectively. We denote r a J to be the sonic point for FGRSF and r\ to be the sonic 
point for flow in ith potential <&i, and define Ar\ as: 

r i _ „gr 

Ar c = (11) 
' c 

In Fig. 2., we plot various Ar*'s (along the z axis) as a function of specific energy £ (along the x axis) and adiabatic 
index 7 (along the y axis) of the flow. The difference in the general profile of variation of Ar l c for various <&i's with £ 
and 7 is quite obvious from the figure. It is observed that when 7 is kept fixed at its minimum value (7 = |), Ar\ 
(Ar c for flows in $1) is always positive (which means r\ > r-F always) and is at a minimum for low energy and starts 
increasing nonlinearly with an increase of £ . At a certain value of £ , it produces a peak and starts falling, with further 
increase of £. The peak value of Ar* is observed to be ~ 0.15. Unlike ArJ, Ar 2 falls very sharply with £ and shows 
a 'dip' at £ ~ 0.04, then starts increasing with £ monotonically and nonlinearly, keeping Ar^ negative for all values 
of £ which indicates that r 2 c < rf r for ultra-relativistic flow (7 = |), □ for all values of £ we consider here. For <I>3 
and $4, we observe that Ar^ and Ar* monotonically decreases with £ and also Arj? and Ar* are always negative e.g., 
r\,r\ < r 9 , r for all £ we consider here. It is observed that for ultra-relativistic flow, $1 is the best approximation to 
study transonic spherically symmetric polytropic accretion to produce the sonic point closest to the sonic point formed 
for FGRSF in general, especially in high energy regime, i.e., for high value of £, while $2 is the worst approximation 
for studying the same phenomena for same set of boundary conditions. If we define a hypothetical quantity Q l rc which 
is a measure of 'goodness' of a particular ith potential $j regarding the closest approximation of sonic points obtained 
for FGRFS, we find that following sequence holds good for ultra-relativistic flow : 

Q\ c > Gtc > Ql c > Q 2 rc (12) 

The situation starts changing as the value of 7 increases. The general profile of Ar\ is unchanged but its peak 
starts shifting towards the lower £ values with increasing 7 and the overall deviation starts increasing. For purely 



nonrelativistic flow (7 = |), Ar 



is found to be quite large. For all values of 7, Ar\ is found to be positive, which 



means that r\ > rf r for all values of 7 and £ . 

However, the behaviour of Ar 2 starts changing drastically with an increase in 7. The 'dip' produced by Ar 2 for low 
values of 7 starts smearing out with an increase in 7 and Ar 2 becomes positive with higher £ for a fixed value of 
7. For higher 7 it has been observed that if we study variation of Ar 2 with £ and 7, Ar 2 starts with a very high 
negative value at lower £ with \Ar 2 \ i max (where |Ar^| maa; stands for the maximum value of the modulus of Ar*), 
monotonically decreasing with 7, reaches zero and becomes positive and increases nonlinearly up to a maximum value 
(this time, a 'peak') and then starts falling again with an increase of £ . So for high values of 7, smearing out of the 'dip' 
in Ar 2 vs £ curve is compensated by the appearance of a 'peak' at some value of £. It is also observed that the peak 
starts shifting towards the higher value of £ as 7 is increased. It is also observed that as we go from ultra-relativistic 
flow towards purely non-relativistic flow, a zero appears in Ar 2 vs £ curve, which indicates that at least one value of 
£ (for a fixed 7) is available where r 2 would be exactly equal to rf r . 

As 7 is increased, not much change is observed in the general profile of Ar\ except that |Arf| decreases with 
increase in 7. The change in profile of Ar\ with increasing 7 is much more interesting. It is observed that with 
increasing 7, not only does a zero appear in Ar\ vs £ (for a fixed value of 7) curve (which means that for some values 
of £ and 7, r\ can be exactly equal to r 9 c r ), but also the value of |Afg| decreases considerably with an increase 
in 7 and as the flow approaches its purely non-relativistic limit, we observe that the overall deviation produced by 
<3?4 in approximating r 9 J is quite small in general (percentage deviation produced by $4 is found to be within the 
limit of 10% for purely non-relativistic flow), which indicates that use of $4 is the best possible approximation of the 
Schwarzschild metric for nonrelativistic flow (flow with high 7 value) , regarding reproduction of the sonic point for a 
flow with fixed £ and 7. 



4-2. The integral curves of motion 

We have mentioned earlier that study of the integral curves, i.e., the variation of Mach number of the flow M with 
radial distance r (measured from the accreting hole in units of r g ) is essential to investigate the 'transonicity' of the 

2 Hereafter, we will describe the flow to be ultra-relativistic for 7 = | and purely non-relativistic for 7 = f according to 
standard practice (Frank, et. al (1992)). 
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loglO(r) 

Fig. 3: Variation of Mach number M (plotted along the Y axis) with radial distance (plotted along X axis in logarithmic 
scale) for accretion and wind for a fixed energy (£ = 0.0123) and adiabatic index (7 = |) of the flow. The solid lines represent 
the variation obtained for FGRSF. For same £ and 7, sonic points r*'s obtained for flow in defferent pseudo-potentials and 
for FGRSF are all different which is clearly observed from the figure. While the value of rj? r is obtained as 32.52 r g , values of 
rl,r^,rl and are obtained as 33.67 r g , 25.27 r g , 31.32 r g and 32.12 r g respectively where r"s are the sonic points for various 
$i's. 



flow. Let us first consider FGRSF. From eqn. 7(c), one can obtain the sonic point rf for a fixed value of £ and 7. The 
value of (^) at rf is then calculated using eqn. 7(d). One can numerically solve eqn. 5(a) and 5(b) simultaneously to 
obtain the value of the Mach number (and other dynamical quantities) as a function of radial distance for a fixed value 
of £ and 7. It is well known that for spherically symmetric accretion onto a Schwarzschild black hole, two solutions 
are obtained while solving the governing conservation equations, e.g., equation for conservation of specific energy and 
from baryon number (or specific entropy) conservation equation. One solution out of these two corresponds to the 
accretion process and the other is for the wind. Using eqn. 6, it is easy to calculate the flow temperature at all points 
of the flow (for accretion as well as for the wind branch). 

For accretion in pseudo potentials $i's, the procedure is exactly the same. Sonic point r l c is obtained by solving eqn. 
10(c) for the ith potential $i for a fixed value of £ and 7. (^) for that particular potential is then obtained by solving 
eqn. 10(d). Starting from r*, eqn. 9(a) and eqn. 9(b) (alternatively, eqn. 9(a) and eqn. 9(c)) could be simultaneously 
solved (using eqn. 10(a)) to get the variation of Mach number of the flow (as well as of other dynamical quantities) 
as a function of radial distance measured in units of r g . The solution for accretion as well as for wind branch can be 
obtained to investigate the 'transonicity' of the flow by plotting Mach number M as a function of radial distance r. 
Temperature of the flow T l (for any ith potential <f>i) can easily be obtained by solving eqn. 6 for di,pi and pi for a 
fixed value of £ and 7. 

In Fig. 3, we plot the variation of Mach number as a function of the radial distance (in logarithmic scale) for FGRSF 
as well as for flows in all $j's. The energy of the flow (for FGRSF and for flows in all $j's) are kept constant at a 
value 0.0123 and value of 7 is taken to be |. It is observed that sonic points of FGRSF r 9 c r , as well as flows in other 
potentials r* are all different, which is normally expected (except for the very special value of £ for which either the 
Ar^ or A?t? vanishes, see Fig. 2(b) and 2(d)). While the value of rf r is obtained as 32.52 r g , values of r\,r%,r^ and r\ 
are obtained as 33.67 r g , 25.27 r g , 31.32 r g and 32.12 r g respectively. It is also to be noted that although the general 
profile of the integral curves are almost same (except for flows in $ 2 where the change in curvature is slightly different 
to that of others, especially for the supersonic accretion branch), the curves do not overlap in general. The amount 
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Fig. 4 (a,b): Variation of the absolute value of the slope of the Mach number M (plotted on the Y axis as ) of accretion 
(Fig. 4(a)) and wind (Fig. 4(b)) with radial distance (plotted along X axis in logarithmic scale) measured from the accretor in 
units of r g . All input parameters are same as that of Fig. 3. Solid line represents FGRSF. The slope in M is found to increase 
towards the accretor non-linearly and monotonically for the accretion branch and the space rate of change of M for <!>i is clearly 
seen to be steepest while that of <E> 4 is observed to be the flattest. For wind, the change is not monotonic; rather it produces a 
peak at the subsonic branch. 
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of deviation (from FGRSF) of integral curves drawn for various $,'s is quite prominent for the supersonic branch of 
accretion and the derivation is observed to be not that significant (except for flows in $2) for corresponding subsonic 
branches of accretion, as well as for the supersonic branch of winds. The deviation maximises with the decrease of 
radial distance from the accretion black hole. This is quite obvious because close to the accretor, the general relativistic 
effects would be much more prominent thus the deviation would be higher. Also it is interesting to note that for the 
same energy £ (and even for a fixed sonic point common to both FGRSF and flows in all <&i's which may be obtained 
by properly tunning the corresponding £'s), the maximum value of the Mach number M max of the flow in different 
$i's are very different. One can observe from the figure that there exists a particular value of r (for a fixed £ and 7) 
where the 'crossing-over' of $4 with FGRSF at the supersonic branch of the accretion is quite distinct, that is, the 
Mach number attained for flows in $4 at a particular point is exactly equal to the Mach number produced by FGRSF. 
Also, the general profile for the variation of Mach number of the flow in $4 is extremely close to the profile obtained 
in FGRSF almost throughout the flow except very close (r ~ 8r g in this case) to the black hole. However, no such 
mutual crossing over is observed in between flows in various $j's. 

If ATfe e and ATth are the relative gain in the kinetic energy and thermal energy of the flow respectively, it is easy to 
show that the change in Mach number AM could be approximated as: 



for spherical accretion in all pseudo potentials discussed here. Also it is easy to understand that the Mach number 
profile as well as the rate of change of Mach number at a particular point should be related to the form of the potential 
$i used to study the flow (and in general to the metric used to describe the spacetime). Also M l max , the maximum 
value of M which could be attended for a particular set of £ and 7 for flows in any zth potential <E>i (for the range 
of r shown here), should be related to the nature of the potential used to describe the flow. It is observed that for 
r > 2r g , M % max anticorrelates with Si (the 'relative stiffness' of the ith potential) for z = 1,3,4 but correlates with <Sj 
for i = 2. 

As change in the Mach number is a result of a mutual tug of war in between change in dynamical velocity and 

polytropic sound speed (alternatively, in between the mechanical and thermal energy of the flow, sec cqn. 13), it 

might be interesting to investigate the variation of the slope of the Mach number (which might be considered as 

the measure of the degree of 'transonicity' of the flow) with radial distance r for two different branches of solutions, 

namely, accretion and wind. In Fig. 4(a) and 4(b) we plot the absolute value of (^r) as a function of loglO(r) for 

accretion (Fig. 4(a)) and wind (Fig. 4(b)) for FGRSF as well as for flows in various <I>iS. For accretion we see that 

J J increases nonlinearly and monotonically as the flow comes closer to the black hole, which indicates that for all 

flows, instantaneous changes in kinetic energy of the flow is always greater than the instantaneous change in thermal 

energy, .i.e., ATke > AT t h for all values of r, which is not the case for winds. For winds (see Fig. 4(b)), change of 

I 1 is not monotonic rather it always shows a peak in the subsonic branch for FGRSF as well as for flows in all $i's. 

The appearance of the peak for flows in pseudo potentials is due to the fact that from the close vicinity of the event 

horizon up to a certain distance r p (r* < r* and r^ r < rf always), the local gain in kinetic energy at any point is 

higher than the local thermal energy gain. So up to r p , ATk e > AT t h but after r p , ATke decreases and starts 

falling as the wind approaches to the sonic point. It is easy to understand that as £ is kept constant throughout the 

flow (for accretion as well as for the wind branch), the bulk motion acceleration of the flow is a continuous process 

throughout the accretion towards the black hole but for winds the acceleration process dominates only up to r p and 

gets the major part of the outward acceleration in the region bounded in r g > r > r l p , though the exact physical 

reason behind this is not clearly understood. However, it should be remembered that the above argument holds only 

for flows in pseudo potentials and not for FGRSF; because for FGRSF, the total energy term can not be decoupled 

2 2 
into various counterparts with individual origin e.g., mechanical energy £k e = \ or thermal energy £ t h — -^bf etc. 

It is clear from the figure that the location of r p (r p for flows in $1) is located furthest away from the black hole and 
ri is located closest to the black hole with the following sequence: 



It is also observed that with increase of 7, the sequence in eqn. (14) is maintained but the location of all r p s moves 
towards the black hole. Thus, the regions from where the wind (in various $i's) is accelerated, moves closer to the hole 




(13) 



r 1 < r 3 < r gr < r 4 < r 



-2 
v 



(14) 



as the flow approaches to its non-relativistic limit. If we define A(^)' to be: 




(15) 
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then one can show that for a fixed value of £ and 7, the following sequence are maintained for various <I>iS as: 



A 



dMY 
~dV) 



< 



max: 



A 



< 



max 



A 



fdMY 

{~d7) 



< 



max 



A 



dM\ 1 

~dV) 



max: 



for accretion where 



^(^r") stands for the maximum value of A(^jf ) . Whereas for wind it is observed that: 

max 

i (} :\ I \ 

A 



fdM 



< 



max 



A 



f dM 

\~~d7 



< 



max 



(dMV 

\-dV) 



< 



max: 



fdM 



max 



Also note that for the accretion branch all A(^)* change sign except A(^) 1 . 




loglO(r) 



Fig. 5: Variation of the temperature deviation (AT's) of the accretion flow with radial distance for ultra relativistic (7 = |) flow 
in various $iS. While the difference in flow temperature in between any T 1 (flow temperature for a particular ith potential 
and T gr (flow temperature for FGRSF) , is plotted along the Y axis, the distance from the central accretor (in logarithmic scale 
in units of r g ) is plotted along X axis. The maximum amount of deviation is obtained close to the accretor, which is expected. 
It is clearly observed from the figure that use of $4 and (except very close to the black hole) $1 are the best approximation of 
the FGRSF flow temperature. Almost the same sort of deviation is found for the wind branch (not shown in the figure). 



4-3. Temperature profile of the flow 

One can easily show that at any point of the flow, the flow temp T l for flow in any particular $iS can be re- written as: 

Ti= lim v f (2 7 -l)(g-^) ) (lg) 



7K [ 2 + M l (7 - 1) 

where M % is the corresponding Mach number for the zth pseudo potential. As it has been stated that M % can be 
calculated at each point of the flow, one can compute T % from knowledge of M l s using eqn. (16). After defining AT 1 
to be: 

rpi rrgr 

AT! /'I 7 s ! 
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where T gr could be computed using eqn. 6, one can plot AT 1 as a function of radial distance (in units of r g ) to 
investigate which <I>i produces the corresponding T l closest to T gr at a particular radial distance and for a fixed value 
of the £ and 7. In Fig. 5. we plot AT 1 vs loglO(r) to demonstrate the dependence of AT 1 on radial distance. It is clear 
from the figure that except for T 1 , all other T l s are always higher than T gr , which indicates that all <I>iS except $1 
produce hotter flow than FGRSF in general. For $1, very close to the hole the flow is cooler compared to FGRSF up to 
a certain distance, after which $1 also produces T 1 higher to T gr . The point at which T 1 becomes equal to T gr , gets 
shifted away from the black hole with increasing 7, keeping £ fixed. Also, it is observed that |AT l | maa; increases with 
increasing 7. From Fig. 5 it is clear that throughout its way from infinity to the event horizon, <&i and $4 produces 
very good approximation of FGRSF regarding calculation of flow temperature. However, after a certain point very 
close to the black hole (after the crossing the point where T 1 = T gr ), $1 produces a sudden and sharp deviation of 
the flow temperature from FGRSF, physical reason for which is not clearly understood. More or less the same sort of 
profile in the AT 1 vs loglO(r) curve also is observed for wind branch. 

5. Conclusion 

In this paper we have solved a set of algebraic and differential equations governing various dynamical and thermo- 
dynamic behaviouars of Bondi (1954) type accretion in a full Schwarzschild metric as well as for motion under a 
number of 'pseudo-Schwarzschild' potentials, to examine the suitability in application of those potentials in investi- 
gating spherically symmetric transonic accretion onto a non- rotating black hole. We have shown that though the 
potentials discussed here were originally proposed to mimic the relativistic effects manifested in the disc accretion, it 
is quite reasonable to use most of the potentials in studying various dynamical as well as thermodynamic quantities for 
spherical accretion. Also, we have shown that depending on the chosen regions of parameter space spanned by specific 
energy £ and adiabatic index 7 of the flow, one potential may be important than others and we could identify which 
potential is the best approximation for FGRSF for what values of £ and 7. We have restricted ourselves to the study 
of simple polytropic flows only. However, the validity of using all these $jS discussed here can easily be examined for 
isothermal accretion and wind as well as for flows with other equations of state. Work in this direction is reported 
elsewhere (Sarkar & Das, 2001). 

It is observed that among all pseudo potentials, <&i (potential proposed by Paczynski and Wiita (1980)) and $4 (one 
of the potentials proposed by Artemova et. al. (1996)) are in general the best in the sense that they provide very rea- 
sonable approximation to the full general relativistic solution. While $1 is the best approximation for ultra-relativistic 
flow, $4 happens to be the best approximation as the flow tends to be fully non relativistic, i.e, 7 tends to have 
the value |. Also we see that there are certain cases for which one or more of the pseudo potentials may give the 
exact match with FGRSF for a particular value of £ or 7 (for a fixed r) in finding some dynamical (r c ,M etc.) or 
thermodynamic (flow temperature T, for example) quantity. 

It is worth mentioning that as long as one is not interested in astrophysical processes extremely close (within 1 — 2 r g ) to 
a black hole horizon, one may safely use the 'pseudo' potentials discussed here to study spherically symmetric accretion 
on to a Schwarzschild black hole with the advantage that use of these potentials would simplify calculations by allowing 
one to use some basic features of flat geometry (additivity of energy or de-coupling of various energy components, i.e., 
thermal (^rj) Kinetic (^) or gravitational ($) etc.) which is not possible for calculations in a purely Schawarzschild 
metric. Also, one can study more complex many body problems such as accretion from an ensemble of companions 
or overall effeciency of accretion onto an ensemble of black holes in a galaxy or for studying numerical hydrodynamic 
or magnctohydrodynamic flows around a black hole etc. as simply as can be done in a Newtonian framework, but 
with far better accuracy. However, one should be careful in using these potentials to study spherically symmetric 
accretion because of the fact that none of the potentials discussed here are 'exact' in a sense that they are not directly 
derivable from the Einstein equations. These potentials could only be used to obtain more accurate correction terms 
over and above the pure Newtonian results and any 'radically' new results obtained using these potentials should be 
cross-checked very carefully with the exact general relativistic theory. 

Although the theory of disc accretion has priority over spherical accretion because of the fact that accretion discs 
describe more realistic situations found in nature, it is not unreasonable to concentrate on spherical accretion because 
for certain cases, that may be quite useful and use of these potentials makes a complicated problem simpler to study. 
For example, for a supcrmassive black hole immersed in intergalactic space in such a way that matter falling on to it 
has negligible intrinsic angular momentum, the accretion (at least close to the hole) is quasi spherical and transonic 
spherical accretion might be a good approximation to mimic the situation. Same sort of approximation is valid when 
an accreting black hole is embedded in a number of donor stars (or star clusters) where the angular momentum of 
the stars arc randomly oriented in such a way that the vector sum of the intrinsic angular momentum carried by the 
accreting matter as a whole may be quite negligible, so as to make Bondi-type accretion a good approximation. In 
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fact, a number of recent works (Coker & Markoff 2001 and references therein, Das, 1999, 2000, 2001a, b, Toropin, et. al, 
1999, Kovalcnko & Ercmin 1998, Titarchuk, et. al. 1996, 1997, Wang & Sutherland 1997, Zampieri, et. al. 1996, Yim 
& Park 1995, Markovic 1995, Tsuribe, et. al. 1995, Kazhdan & Murzina 1994, Fortncr 1993) still deal with spherical 
accretion to investigate some basic astrophysical processes the black holes and neutron stars. So we believe that work 
presented in this paper is relevant and will be useful in investigation of various aspects of accretion and wind around 
non-rotating and slowly rotating compact objects. 
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